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ABSTRACT 


The NAVSTAR Global Positioning System (GPS) has provided a quantum leap in 
real-time autonomous navigation capabilities. NASA’s Space Shuttle will be receiving an 
integrated GPS capability in the near future, and the orbiter Endeavour has been equipped 
with a stand-alone GPS receiver. Although much data is available regarding spacecraft 
GPS receiver performance at higher altitudes, little information is available for spacecraft 
at Shuttle altitudes of approximately 400 km where drag and gravity effects are more 
pronounced. 

GPS receiver navigation solution data from Shuttle mission STS-69 was made 
available by NASA and provided an opportunity for evaluating GPS performance in low 
Earth orbit. NASA ground tracking network and Tracking and Data Relay Satellite 
(TDRS) data for this mission provided a reference for comparison. Analysis of the data 
was accomplished using Satellite Tool Kit (STK) for visualization and Matlab routines for 
data comparison. 

GPS navigation solutions were available for approximately 65 percent of the STS- 
69 mission, and they generally coincided with the reference track. Differences between the 
GPS navigation solution state vectors obtained using the Standard Positioning Service 
(SPS) and the reference state vectors produced RMS position differences between the 
data sets of about 1500 m. One sigma position accuracies of 54 m in the vertical direction 
and approximately 1400 m in the downtrack direction were experienced. Velocity vector 
magnitude differences during this period were generally + Im/s, with a RMS velocity 
difference of less than 9 m/s. One sigma velocity accuracies of approximately 4.2 m/s in 
the vertical direction, 2.3 m/s in the downtrack direction and 1.5 m/s in the crosstrack 
direction were experienced. A firm conclusion regarding Shuttle GPS accuracies could 
not be drawn because all sources of error were not identified. Based on these results GPS 
appears to be an excellent navigation source for Shuttle state vector information; however, 
another navigation source such as INS must be present to provide a check against 


spurious data points and periods of outage. 





TABLE OF CONTENTS 


ДУ 1Ш1010К011 21 ааа... 1 
ИОК Г К _...........................1Ш.................................. 1 
TT ODE LOI 1 

СМ ЫН тос ТЭ? 0000000000 n B 

Е ЛОВИТИ п. 010.00 аа ааа... 6 
РЕ iere teretes robe te eee aeta nnne nnn ntn nen 8 
П. REFERENCE FRAMES AND SHUTTLE TRACKING .............1. 1... 11 
AAN SA ——o————-——«— 11 
| еТкетпа! Кетегепсе Етатев. .......................................................... 11 

2. то Г ТОЭЦЭНУЭ 0 Цэ 5277 оэотэ тээг анална, a 

3. Transformation {о Ј2000 Кеѓегепсе Егате ................................................ 13 

Пр ао епенсе гате 0... ааа ааай 13 

a. ЫГ ЭЭЭ ЭТГ9777...... WB enc tti tes o ЙГ. 15 

Eo SN ЕРИ АСК МА аанак. 16 
IAE TES TOOLS AND PROCEDURE а. 17 
ИКЕ MA ШЕК КККК .................................................................... 17 
т Махаапоп И В... дА... 17 

Lo IAE 2 .................................................. 17 

САКЕ е CCE 2 ЕСО... аана. 18 
о. ПИ И ПВ и 18 
O nia nadaa н... 18 
Е dm nda en E 18 
Шошан ы аы наанаа. Do 

A. A аео саа 2? 

ИМ е 7 277 TOLL ECCE 22 
И. а... ое. 23 

RR 22 22222 еее 25 
ога ага Piles О п... ағас... + 23 

2. Edit GPS Files with Unix Scripts and Ёоптгап Ргортат.............................. 23 
ВОЕН ОВ WILD ЛИК ЭС... 23 

ЕИМШ ен сг ги ОШГап Ргботап! _...................................................... 26 

c. Reference Frame Transformation Using ConvertTool .................... 27229 

БМ Бор 0 222 ............................ 29 

ШЕШ РО расаеа Паја Ке "пи ааа. 31 

ЕМН ӨГ 110) guae m T TT I T TT ae: 31 
агланав еее, 31 
ЕНЕР 22222272 анаа. 33 
ЕМИ ШЕ ОЕЕО БОМ КЕЗО oe aaa ieaiaia 33 
О НИНИН ЕМ. 33 
DA 0 nn. nennen rap 35 
ОЕЕО аана ое ens 36 
Нэгэ. 20:12 2 ___...................................................................... 36 


B. STK RESULTS... нер 58 


C. SHUTTLE UEN COMPARISON RESULTS. nn. 79 

1. ЏЕМ Тгапејогтаџоп. ._._____. мм 79 

2. UEN Results ВОИНИ 80 

V. CONCLUSIONS AND RECOMMENDA TON S ns 95 
A. CONCLUSIONS SS. 95 

В. ВЕСОММЕМРАТІОМЅ ЭЭГ sS 96 
1165Т ОЁ ЁЕЕЕАЕХСБЕ5...... л 22222222 22. 97 
APPENDIX A. FORTRAN PROGRAM FOR EDITING PROPAGATED FILE............ 99 
APPENDIX B. FORTRAN PROGRAM FOR EDITING NAVG-11 STATE VECTOR 
ЕШЕ Т... И 103 
APPENDIX С. STK VEHICLE FILE о. 107 
APPENDIX D. MATLAB M-FILE FOR COMPARING GPS AND PROPAGATED 
DATA 2 0 7 оо 113 
APPENDIX E. MATLAB M-FILE FOR PLOTTING GPS ORBIT IN 3-D............... 119 
APPENDIX F. MATLAB M-FILE FOR COMPARING GPS AND NAVG-11 STATE 
VECTORS WITH MATCHING TIMES eee eei 123 
APPENDIX G. MATLAB M-FILE FOR TRANSFORMING MATCHING GPS AND 
NAVG-11 STATE VECTORS TO UEN FRAME FOR COMPARISON .................. 133 
APPENDIX H. DAY 250 ТКР1ОТ5........ SS 141 
APPENDIX I. РАҮ250 МАТГАВ РГОТЅ 77/5 145 
APPENDIX J. DAY 251 STK PLOTS ............. ssec И 151 
APPENDIX K. DAY251 MATLAB PLOTS ............. ees 157 
APPENDIX L. DAY 252 STK PLOTS „= № 187 
APPENDIX M. DAY 252 MATLAB PLOTS |... ees 193 
APPENDIX N. DAY 253 STK PLOTS ——— и 221 
АРРЕМПІХ О. БАҮ 253 МАТІ.АВРІОТ6...22..2.....4Д4 4ӘО: 2. 231 
APPENDIXP. DAY 254 STK PLOTS 7 аа” и 237 
APPENDIX Q. DAY 254 МАТГАВ РГОТ5. ла. 243 


уш 


ДЕРЕ эг РАУ 255 МАТГАВРЕОТӘ........................................................... 25] 
A IIA ЗС SIN а О 7_.............................................................. 263 
SEBENDISTIZDBANZSSMATEABPLOTS zen... anne. 269 
RENDIDA S AST PETS... erreurs 215 
шингээх БАҮ 257 МАТРБАВРІЕОТЗ.......................................................... 283 
Ат 2-5 ЭТК РЕО)ТУ...................................................................... 289 
ЕВ ЕМО ЮА 255 МАТГАВРЕО15............................................................. 295 
СООО ОЗИ DAY 250 У ПС РЕОГЗУ ла... nn 301 
ИНЕМ АС ПАР 2595 МАТГАВРЕОТЗ..................................2т................ 307 
ВВ АВ. ПРАХ 260 5ІКРІЕОТ5...2-57........................................................ 313 
SEBENDIX AC рАҮ 2ООМАТГАВРІОГЫ .......................................................... 22! 
ПЕЧЕ ШАВИ ОМЕТ ан 221 





ACKNOWLEDGEMENT 


The completion of this effort would not have been possible without the support 
and assistance of several generous people. Sincere thanks is extended to the following: 
My wife Lynn for her untiring support; Professor Sandi Scrivener for her encouragement 
and guidance; Professor Don Danielson and Professor John Clynch for the insights they 
provided; Mike Gabor and Scott Wallace at the Air Force Phillips Lab for all the time 
they graciously spent teaching me about orbit determination and how to manipulate the 
data; Dr. Jim Woodburn of Analytical Graphics for his STK and programming assistance 
and expertise; Ray Nuss and Russell Carpenter from NASA JSC and Rosanne Nikolaidis 
of RSOC for answering an unending stream of questions and providing the crucial data; 
and COL John Casper, CDR Kent Rominger, Dr. Andrew Thomas, and MAJ Michael 
Knapp of the Astronaut Office for their ideas, guidance, leads and points of contact. 





E 


I. INTRODUCTION 


The NAVSTAR Global Positioning System (GPS) has provided a quantum leap in 
real-time autonomous navigation capabilities. GPS technology has been applied in 
numerous fields and has now become a method for autonomous spacecraft navigation. 
Flight experience with the Topex/Poseidon spacecraft has indicated that very precise 
spacecraft orbit information can be obtained based on GPS measurements. NASA’s 
Space Shuttle will be receiving an integrated GPS capability in the near future, and the 
shuttle Endeavour has been equipped with an operational stand-alone system. Although 
much data is available regarding spacecraft GPS receiver performance at higher altitudes, 
little information is available for spacecraft at Shuttle altitudes of approximately 400 km 
where drag and gravity effects are more pronounced. 

GPS receiver navigation solution data from Shuttle mission STS-69 has been made 
available by NASA and provides an opportunity for evaluating GPS performance in low 
earth orbit. NASA ground tracking network and Tracking and Data Relay Satellite 
(TDRS) data for this mission was also available for comparison. The objective of this 
thesis was to compare Shuttle GPS receiver navigation solutions with the NASA reference 
track. Data comparison routines written in Matlab were used to analyze both sets of state 
vectors. The visualization features of the orbit analysis software package Satellite Tool 


Kit (STK) aided the analysis of the data sets. 


A. GPS BASICS 


1. System Operation 


GPS was developed by the Department of Defense (DoD) as a worldwide satellite- 
based radionavigation system which provides users with position, velocity and time 
information. Using the concept of ranging, the system measures the distance to several 
satellites to determine position. Ranging is dependent on time synchronicity among the 


receiver clock and the satellites. The satellites transmit radio pulses in the form of distinct 
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pseudo-random sequences at specific known times. The receiving equipment recognizes 
these sequences and measures the precise time at which the pulses arrive. Based on the 
transmission and reception times, the time it took the pulses to travel from the satellite to 
the receiver can be determined. Since the speed of light is constant and the travel time of 
the pulses is now known, the range to the satellite can be calculated. 

Determining the range to a particular satellite places the user on a sphere with a 
radius equal to this range. The intersection of two spheres of position is a circle, and the 
intersection of three spheres is two points. Position can be determined using three satellite 
ranges because the ambiguity of the second point can usually be resolved as an 
unreasonable solution. The intersection of four spheres of position produces a single 
point. (Herning, 1996, p. 46) 

Ranging requires that the receiver’s clock be synchronized with the satellite’s 
clock; however, in order to support receivers with less than perfect chronometers, it is 
necessary to apply correction and this requires a fourth satellite as shown in Figure 1.1. 
First, ranges to the four satellites are calculated using the original receiver clock time and 
are called pseudo-ranges. The resultant four spheres of position would intersect in a 
single point if there were no clock error. The receiver clock error is then determined 
algebraically. Because there is only one value of clock error for which the four spheres 
can be made to intersect at a single point, the user’s position solution is thus provided as 
shown in Figure 1.2. (Herring, 1996, pp. 46-47) 

Arranged in six orbital planes with four satellites in each plane, the constellation 
of GPS satellites consists of twenty-four satellites. This ensures that a minimum of four 
satellites are always observable by a user anywhere in the world. The satellites are in 
circular orbits inclined at fifty-five degrees with altitudes of 20,200 km and complete an 


orbit every twelve hours. 


aw UNCORRECTED SPHERE 
OF POSITION 


CORRECTED SPHERE 





Figure 1.1. Clock error in the receiver causes the range to the GPS satellite to be 
incorrect. The resulting spheres of position (thick lines) do not intersect at a single point. 
By adjusting the receiver clock the ranges are corrected and meet at one point (thin lines). 
In this figure this appears to require three satellites; however, in three dimensions four 
satellites are required. From Herring, 1996, p. 46. 
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Figure 1.2. Navigation Using GPS. From NAVSTAR-GPS Joint Program Office, 1987, 
p. 1-17. 


2 GPS Accuracy 


GPS provides two levels of service, the Standard Positioning Service (SPS) and 
the Precise Positioning Service (PPS). Utilizing spread spectrum techniques, GPS 
satellites transmit on two L-band frequencies, L1 at 1575.42 MHz and L2 at 1227.6 MHz. 
Two spreading functions known as C/A-code, or Coarse/Acquisition code, and P-code, or 
Precise code, are employed. Currently only the C/A-code signal is available to all GPS 
users and provides the SPS, while the P-code is available to DoD-designated users only 
and is necessary for the PPS. Although user position can be determined using both codes, 
only the P-code is transmitted on both frequencies. This enables users with P-code access 
to perform a dual frequency comparison to compensate for ionospheric effects while C/A 
code users must rely on a less accurate model of the ionosphere. Satellite ephemeris 
(orbital element) data, almanac data for finding other satellites in the constellation, 
atmospheric propagation correction data, satellite clock-bias information, system time and 
satellite status information are provided by the navigation message which is transmitted on 
both frequencies. (U.S. Naval Observatory, 1996, pp. 1-2) 

Full accuracy is denied to non-DoD designated users through the use of Selective 
Availability (SA), which modifies the navigation message. SA reduces accuracy by 
altering the satellite ephemeris data and clock frequency, a procedure which is known as 
dithering. Anti-spoofing protects against false satellite transmissions by encrypting the P- 
code, forming the Y-code. (U.S. Naval Observatory, 1996, p. 2) 

The SPS has a 95 percent probability of providing horizontal positioning accuracy 
within 100 meters, vertical positioning accuracy within 156 meters and timing accuracy 
within 340 nanoseconds. (U. S. Naval Observatory, 1996, p. 1) Typically, one sigma 
positioning accuracy is on the order of 50 meters horizontally and 75 meters vertically 
with a velocity accuracy of 50 cm/s. The specification for PPS requires positioning 
accuracy of 16 meters Spherical Error Probable (SEP) in contrast to the SPS specification 
of 76 meters SEP and a velocity accuracy of 10 cm/s. A positioning accuracy of 10 
meters SEP is usually experienced. (Clynch, 1996) Sources of GPS inaccuracy and their 


impact are shown in Table 1.1. 


Component 
Receiver Clock 1000 km 
Selective Availability 30m 


Ionosphere 1-30m 


Atmosphere 1-6 тш 
Orbit Error 1-3 m 
Satellite Clock 1-3 m 
Multipath - C/A Code 0.5-150 m 
Multipath - P Code 0.1-15 m 


Receiver Noise 0.1-1m 





Table 1.1. GPS Sources of Inaccuracy. From Clynch, 1996. 


GPS system time is maintained by the Composite Clock which incorporates all 
operational Monitor Station and satellite frequency standards. This system time is 
referenced to the U.S. Naval Observatory’s Master Clock within a standard tolerance of 
one microsecond. Over the last several years GPS system time has maintained a tolerance 


within a few hundred nanoseconds of the Master Clock. (U.S. Naval Observatory, 1996, 
p. 3) 


B. SHUTTLE GPS 


An integrated GPS navigation capability has been developed for the Shuttle. The 
scheduled replacement of TACAN air navigation stations has been the primary motivating 
factor in developing this capability. Currently TACAN provides the primary Shuttle entry 
navigation aid from post blackout to final approach. During final approach through 
landing the Microwave Scanning Beam Landing System (MSBLS) is the primary 
navigation aid. Presently, only the orbiter Endeavour is equipped with GPS equipment. 
The other orbiters will be ёоо as they are upgraded. 

Two options were considered for integrating GPS into the Shuttle navigation 


system. The first option was to augment Shuttle entry navigation filters to process GPS 


measurements. The second option was for the Shuttle navigation system to process state 
vector outputs from the GPS receiver. Either option can be extended to other Shuttle 
mission phases such as ascent and on-orbit operations. The second option of processing 
GPS state vector outputs was selected because it required fewer software changes and had 
the minimum impact on current system operation. Additionally, it facilitated the design 
objective of a phased-development approach which would permit evaluation of GPS 
receiver performance during flight before it was used as the primary entry navigation aid. 
(Kachmar et. al., 1993, pp. 313-317) With respect to developmental testing for STS-69, 
the GPS receiver navigation solution was not available to the orbiter in real time even as a 
back up. Beginning with STS-79, the navigation data will be processed by the orbiter 
computers and real-time state vectors will be available during launch and landing. 

Using the second approach, a GPS receiver state vector can be used to reset the 
baseline navigation state. This navigation state reset procedure is similar to a ground state 
vector uplink with the important difference that the GPS state vector quality assessment is 
done onboard the Shuttle. In order to be selected for guidance, the GPS state vector must 
meet a predefined transient tolerance level. Transients can be induced by the Shuttle 
mission phase environment, the GPS receiver switching between satellites, and periods of 
satellite reacquisition. The GPS receiver can be provided aiding data by the Shuttle's 
Inertial Measurement Units (IMU) to assist the receiver in determining its position for 
acquisition of the GPS satellites. (Kachmar et. al., 1993, p. 317) 

The Shuttle was designated as a DoD-authorized user of PPS due to the need to 
guarantee a minimum level of entry navigation performance during times of national 
emergency. It was also desired that current off-the-shelf receiver designs be used that 
complied with the GPS Joint Program Office Requirements Document. (Kachmar et. al., 
1993, p. 317) The Collins 3M receiver was selected for installation onboard Endeavour. 
The five channel 3M receiver tracks four satellites simultaneously while reading the 
navigation message for the next satellite to be used by the receiver. It is a P-code 
receiver with position accuracy of 15 meters SEP, velocity accuracy of 10 cm/s, and 
timing accuracy of 40 nanoseconds. (Collins, 1994) (Nuss, e-mail, 1996) Two antennas 


are employed to provide maximum coverage for the receiver. The upper-hemi and lower- 


hemi antennas are located on the top and bottom of the orbiter respectively as shown in 


Figure 1.3. 
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Figure 1.3. Shuttle GPS Antenna Placement. From Carpenter and Hinkel, 1995, p. 6. 


С: STS-69 


Shuttle mission STS-69 lifted off from Pad 39-A at the Kennedy Space Center 
(KSC) on September 7, 1995, at 11:09:00.052 am. EDT. Endeavour’s crew of five 
included: Commander David M. Walker, Pilot Kenneth D. Cockrell, Payload Commander 
James S. Voss, Mission Specialist James H. Newman, and Mission Specialist Michael L. 
Gernhardt. The mission featured the second flight of the Wake Shield Facility (WSF), a 


four-meter-diameter saucer-shaped free-flying satellite designed to grow semiconductor 
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films in the ultra-vacuum created in the wake of the satellite as it moves through space. 
The Spartan 201 free-flying astronomy satellite was also deployed and retrieved. 
Additionally, the crew tested assembly techniques for the international space station and 
tested thermal improvements made to space suits used during space walks. On September 
18, 1995 at 7:37:56 am EDT, after 10 days, 20 hours, 28 minutes and 55 seconds, STS-69 
landed at KSC. After a successful mission and traveling over four and a half million miles, 
Endeavour touched down on Runway 33 at the KSC Shuttle Landing Facility. 

The WSF was employed to conduct several GPS experiments sponsored by the 
Texas Space Grant Consortium and was equipped with a TurboRogue GPS receiver 
furnished by the University Corporation for Atmospheric Research (UCAR). University 
of Texas at Austin researchers and Johnson Space Center personnel conducted the 
experiments which were aimed at: assessing relative positioning using the TurboRogue 
receiver on the WSF and Endeavour’s Collins 3M receiver, evaluating precision orbit 
determination in a low-altitude environment (400 km) where perturbations due to 
atmospheric drag and the Earth’s gravity field are more pronounced than for higher 
altitude satellites such as TOPEX/POSEIDON, and determining atmospheric temperature 
profiles using GPS signals passing through the atmosphere. (Schutz et. al., 1995, pp 1-2) 


During STS-69 Endeavour’s 3M receiver operated in a SPS mode. 





П. REFERENCE FRAMES AND SHUTTLE TRACKING 


A. REFERENCE FRAMES 


Several reference frames were used when comparing GPS navigation solutions 
with the reference track. Shuttle navigation and tracking sources utilize the Arıes Mean- 
of-1950 Earth-Centered Inertial reference frame, Greenwich Mean Time and English units. 
GPS navigation solutions provide position data ın the WGS 84 Earth-Centered Earth- 
Fixed reference frame using Coordinated Universal Time and metric units. GPS velocity 


data is in a Shuttle Up-East-North frame and in metric untis. 


1. Inertial Reference Frames 


The Aries Mean-of-1950 (M50) and J2000 Earth-Centered Inertial (ECI) reference 
frames are not actually fixed with respect to the mean position of the stars in the vicinity 
of the Sun, or inertial space. These inertial frames are defined by the direction of the 
Earth’s axis of rotation, or celestial pole, and the direction from the Earth to the Sun on 
the first day of spring which is known as the vernal equinox or first point of Aries. The 
celestial pole is affected by gravitational forces, primarily those exerted by the Sun and 
Moon. These forces produce a quasi-conical motion of the mean celestial pole around the 
pole of the ecliptic, or axis of rotation of the Earth’s orbit around the Sun. This lunisolar 
precessional motion has a period of 26,000 years. The motion of the actual celestial pole 
around the mean pole is known as nutation and has a period of 18.6 years. (Larson and 
Wertz, 1992, pp. 94-95) (Seidelmann, 1992, p. 12) The International Astronomical Union 
(IAU) has adopted conventions for calculating general precession (IAU 1976) 
(Seidelmann, 1992, p. 103) and general nutation (TAU 1980) (Seidelmann, 1992, p. 111); 
however, they do not account for transient pole wander. 

Motion of the vernal equinox subsequently accompanies the motion of the celestial 
pole. As a result, in order to define an ECI reference frame the date which specifies the 
position of the vernal equinox must be established. (Seidelmann, 1992, p. 12) The J2000 


frame is a right-handed Cartesian coordinate system with its origin at the center of the 


1| 


Earth and its X and Z axes pointing towards the mean vernal equinox and mean rotation 
axis of January 1, 2000. The M50 coordinate system uses the mean vernal equinox and 
mean rotation axis of 1950 as its reference directions. A constant transformation matrix 
has been developed to perform the conversion from M50 to J2000 coordinates. 


(Seidelmann, 1992, pp. 184-187) 


2 WGS 84 Reference Frame 


The World Geodetic System 1984 (WGS 84) reference frame is a Conventional 
Terrestrial System (CTS) whose ongin is the center of mass of the Earth. The system’s Z- 
axis is parallel to the direction of the Conventional Terrestrial Pole (CTP) defined by the 
Bureau International de l'Heure (BIH). The X-axis is the intersection of the WGS 84 
reference meridian plane, which is parallel to the Zero Meridian defined by the BIH, and 
the CTP”s equatorial plane. The Y-axis completes a right-handed, earth-fixed orthogonal 
coordinate system. This reference frame rotates with the Earth and assumes that the earth 
is rotating at a constant rate around a mean celestial pole, the CTP. Variations from this 
approximation cause the WGS 84 CTS to differ from the true Instantaneous Terrestrial 
System (ITS) which is rotating around an instantaneous pole, the Celestial Ephemeris Pole 
(CEP). The WGS 84 CTS can be related to the J2000 Conventional Inertial System (CIS) 
by Equation 2.1. 


CTS - [A] [B] [C] [D] CIS (2.1) 


The rotation matrices account for: [A] polar motion, [B] sidereal time, [C] 
astronomic nutation, and [D] precession. Matrices B, C and D establish the relationship 
between the CIS and the ITS. The ITS is related to the CTS by matrix A which provides 
the relationship between the CEP and the CTP. (DMA, 1987, pp. 2-1 - 2-3) 


12 


3: Transformation to J2000 


The Shuttle navigation system and the NASA Shuttle tracking sources produce 
state vectors using the M50 ECI reference frame, and GPS state vectors are provided in 
the WGS 84 Earth-Centered Earth-Fixed (ECEF) reference frame. As shown in Figure 
2.1, the MSO frame is fixed in inertial space while the WGS 84 frame rotates with the 
Earth. In order to compare state vectors from both systems it is necessary to perform a 
transformation between the two reference frames. This transformation requires an 
intermediate step of converting the data into the J2000 ECI reference frame. For this 
reason comparison of state vectors from the two sources was conducted in the J2000 


frame. 





Vernal Equinox 
Direction 





Celestial Coordinates Earth-fixed Coordinates 


Figure 2.1. ECI and ECEF Coordinate Systems. From Larson and Wertz, 1992, p. 95. 


4. Up-East-North Reference Frame 


GPS velocity data is provided in the Up-East-North (UEN) reference frame which 
has its origin fixed with the spacecraft as shown in Figure 2.2. The ЏЕМ system can be 
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visualized using the idea of the Shuttle flying with its belly facing the Earth much like the 
attitude an aircraft would assume while flying above the surface of the Earth. The 
Shuttle’s nose would point along the velocity vector. Up would be equivalent to altitude 
above the ground, East would point straight ahead from the nose and North would point 
in the direction of the left wing tip. Technically, the Up axis is in the radial direction from 
the center of the Earth. The East axis is in the local orbital downtrack direction and 
parallels the velocity vector. The North axis is in the crosstrack direction completing a 


right-handed coordinate system and parallels the orbit normal. 


Z 
DOWNTRACK 
CROSSTRACK 
U 
2 

QR | 
12% | 
| 
| 

Ү 


Figure 2.2. Up-East-North Reference Frame. 


In order to perform the transformation from the UEN reference frame to the WGS 
84 reference frame two Euler angle rotations must be performed. A rotation about the Y- 
axis by +ò followed by a rotation about the Z-axis by a value of -y is required as shown in 


Equations 2.2, 2.3 and 2.4. 
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X ОД 


у| -С4-у)С/9)|Е (2.2) 
WGS84 N 
X cos(-v) sin(-y) O||cos(ö) 0 -sin(ö)| |U 
Ү =|-sin(-7) cos(-y) 0 0 1 0 Е (2.3) 
Mm 0 0 11! 51(0) 0 cos(ö) | | N 
X cos(y)cos(ö) -sin(y) -сов(у)віп(6)| 10 
ү =|sin(y)cosíó) cos(y) -sin(y)sin(ö) | | E (2.4) 
2 sin(d) 0 со5(9) М 


WGS84 
This transformation is easily realizable in Fortran and was applied to GPS state 


vector velocity components. 


5, Time 


International Atomic Time (TAI) compares several processes of physics such as 
cesium frequency standards and hydrogen masers to form a standard timescale that can be 
used to uniquely identify the instance of time at which an event occurs on Earth. 
(Seidelmann, 1992, p. 2) Universal Time (UT) is based upon the Earth’s rotation with 
respect to the Sun. Since the rotation of the Earth is affected by irregular forces, UT is 
irregular with respect to TAI. The difference between TAI and UT is increasing irregularly 
by approximately one second every eighteen months, and this difference is always 
maintained as an integral number of seconds known as leap seconds. Coordinated 
Universal Time (UTC) is an atomic timescale that is kept in close agreement with UT 
through the addition of these leap seconds. (Seidelmann, 1992, pp. 6-7) UTC is the 
timescale used by GPS and is provided in the receiver navigation solution. | 

Greenwich Mean Time (GMT) is the basis of civil time for the United Kingdom 
and is subsequently related to UTC; however, in navigation terminology GMT has been 
used as UT. Asa result, the term GMT is somewhat ambiguous. (Seidelmann, 1992, p. 7) 
All Shuttle tracking data uses GMT as its timescale. In this specific case, GMT is 


equivalent to UTC. (Nikolaidis, e-mail, June 1996) In order to facilitate the use of the 
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coordinate frame transformation program and STK, all times were converted to Mission 


Elapsed Time (MET) in seconds which starts at the time of launch. 


B. SHUTTLE TRACKING 


NASA uses two primary sources for Shuttle tracking data, ground sites using C- 
band radar and Tracking and Data Relay Satellites (TDRS). Several of the NASA 
tracking sources are listed in Table 2.1. C-band radar sites use range and angle data and 
TDRS uses doppler measurements to create Shuttle state vectors. Normally two to four 
state vectors are generated per orbit, although, more vectors are generated during periods 
of high activity such as rendezvous, translational maneuvers and IMU alignments. With 
the exclusion of rendezvous, the three sigma position accuracy of the state vectors is 200 
m in the radial direction, 450 m in the downtrack direction and 200 m in the crosstrack 
direction. A three sigma velocity accuracy of 0.45 m/s in the radial direction, 0.20 m/s in 
the downtrack direction and 0.25 m/s in the crosstrack direction is typically achieved. 
During quiescent periods, three sigma accuracies of 100 m and 0.30 m/s radially, 250 m 
and 0.10 m/s downtrack and 100m and 0.15 m/s crosstrack may be obtained. During 
rendezvous and other highly active periods errors on the order of several kilometers may 


be experienced. (RSOC, 1996, p 1) 


Source Identifier 
East TDRS 

West TDRS 

Ascension Island 


Bermuda 


Kwajalein Island 
Wallops Island 
Table 2.1. NASA Shuttle Tracking Sources. 





16 


Ш. DATA FILES, TOOLS AND PROCEDURE 


A. DATA FILES 


1. GPS Navigation Solutions 


GPS data files from the 3M receiver for each flight day of STS-69, Julian days 250 
through 260, were made available by the Johnson Space Center (JSC). There was an 
average of 60 navigation solution files per day and a total of 660 files. Each file had from 
800 to 3300 state vectors with fixes obtained generally every second, although every day 
had some period of outage with no GPS data. Individual state vector entries contained: 
GPS system time; Coordinated Universal Time (UTC); position in XYZ WGS 84 ECEF 
coordinates (m); velocity in UEN coordinates (m/s); position in latitude, longitude and 
altitude; pitch, roll and heading information; acceleration in UEN coordinates; and receiver 


health and status information. A typical GPS data file is shown in Figure 3.1. 


2: NAVG-11 State Vectors 


Rockwell Space Operations Company (RSOC) provided a NAVG-11 Shuttle 
Navigation Postflight Product which included ground navigation history vectors. The 
NAVG-11 product contains the orbiter solution vectors which are processed real-time 
using C-band range and angle data and TDRS doppler data. Each entry contained the 
batch ID, orbit number, batch number, UTC labeled as Greenwich Mean Time (GMT), 
and position and velocity in XYZ M50 ECI coordinates (ft and ft/s). There were 
approximately 430 state vectors for the entire mission in this file. Two to four vectors 
were provided per orbit with the bulk of the data coming from TDRS. (RSOC, 1996, p. 
1) A portion of the NAVG-11 file is shown in Figure 3.2. 
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3; Propagated State Vectors 


RSOC also provided a second more densely propagated data file with over 16,000 
state vectors for the entire mission. The first column of the propagated file is Mission 
Elapsed Time (MET) in hours from launch which was given to be 1995:250:15:09:00.0 
GMT. The program that produces the denser ephemeris propagates from vector to vector 
using a Nystrom-Lear fourth order integrator for both position and velocity. Four 
derivative evaluations are required for each 60 second integration step. The program 
takes into account modeled translational maneuvers, Sun and Moon induced gravitational 
perturbations, drag computed using a Jacchia 1970 atmosphere model, and a 7x7 
gravitational potential model. Figure 3.3 shows a small portion of the propagated data 
file. (Nikolaidis, e-mail, March 1996) 


B. TOOLS 


1. File Transfer Protocol 


All analysis was performed on a Silicon Graphics Indigo П workstation. File 
Transfer Protocol (FTP) provided the only manageable means for retrieving the GPS data 
files from the JSC server. The large number and size of the files made the option of 


copying them to some form of media for transfer unattractive. 


2. Unix Scripts 


Unix scripts, sequences of Unix commands that can be stored in files, were 
employed in editing all the data files and were crucial ın the editing of the GPS data files. 
Hand editing all the data files was not feasible, and the Unix scripts proved to be the only 
way to deal with the large quantity of data files by automating the procedure of 
uncompressing and compressing files, grouping files by day, editing files and executing 


Fortran programs. More specifically, the Stream Editor (SED) command made the 
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It provided an automated means for 


analysis of the enormous amount of data possible. 


headers and troublesome characters such as colons and slashes. 


Ф 
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Figure 3.1. GPS Navigation Solution Data File. 
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Figure 3.2. NAVG-11 State Vector Data File. 
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Figure 3.3. Propagated Data File. 
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3. Fortran Programs 


Fortran programs proved to be very versatile and effective in dealing with the 
varying data files. They were used to select the desired state vector elements from the 
input file line entries, convert the input file times to MET in seconds, and write the new 
data files. In the case of the GPS data files, Fortran programs were also used to thin the 
dense GPS data files by a factor of ten and transform GPS UEN velocity data into WGS 
84 coordinates. Fortran routines also converted all the NAVG-11 and propagated data 


from English units to metric units. 


4. Satellite Tool Kit 


Satellite Tool Kit (STK) is an orbit analysis software package developed by 
Analytical Graphics that employs an object paradigm. A Graphical User Interface (GUI) 
is used to manipulate objects such as scenarios, vehicles, facilities, targets, and sensors. 
This thesis required work with only scenarios and vehicles. A vehicle is defined as a 
movable land, sea, air or space object. A scenario is a set of objects that is to be 
visualized. The objects in the scenario are related by time to a map background which 
provides a geographic reference frame. (Analytical Graphics, 1996, pp. 1-2) STK Vehicle 
files for the Shuttle were created using GPS ephemeris data and the propagated file 


ephemeris data and incorporated into scenarios for visualization. 


5, ConvertTool 


The STK math utility convertTool provided a means for transforming data files 
from one reference frame to another. In particular it permitted files to be transformed 
from M50 to J2000, J2000 to M50, ECEF to J2000, and J2000 to ECEF. The program 
accounts for general precession according to [AU 1976 and general nutation according to 
[AU 1980. It does not account for the Earth’s irregular rotation rate or pole wander. 


(Woodburn, e-mail, 1996) 
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6. Matlab 


Matlab is an interactive numerical computation, data analyis and plotting software 
package. The basic Matlab data element is a matrix, and it uses matrix manipulation to 
perform many numeric calculations in a fraction of the time it would take to write and run 
a program in a language such as Fortran which would perform the same function. Matlab 
allows the execution of sequences of commands that can be stored in files called M-files. 


Matlab was used extensively to compare GPS data and Shuttle reference track data. 


€. PROCEDURE 


Upon undertaking a project with such a large amount of data, procedures needed 
to be established to facilitate the transfer and manipulation of the vast data sets. From the 
beginning, using File Transfer Protocol to retrieve the data, to the end, using Matlab to 
perform major coordinate transformations, the emphasis was on using computer programs 
and tools to simplify all aspects of working with the data. The procedure and tools 


utilized in working with the data is outlined in the flow chart in Figure 3.4. 


1. Obtain Data Files 


File Transfer Protocol (FTP) was employed to retrieve the 660 GPS data files from 
the JSC server. The more manageable NAVG-11 and propagated files were obtained on 
floppy disk and 8 mm tape, respectively. Grouping the files by flight day provided the 
most efficient way of handling the data. The best results were achieved by plotting each 


day’s worth of data in both STK and Matlab. 


2: Edit GPS Files with Unix Scripts and Fortran Program 


a. Automation with Unix Scripts 
The next step was to edit and manipulate the files using Unix scripts and 


Fortran programs. The following Unix script, “thinfiles,” was invoked to execute the 
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scripts and the Fortran program that would manipulate all 660 GPS data files for flight 
days 250 through 260. 
thinfiles 


fileOscript 
filelscript 
file2script 
file3script 
file4script 
fileSscript 
file6script 
file7script 
file8script 
file9script 
file60script 


Each file#script in “thinfiles” executed the “edit1script” script for all the GPS data files in 


each flight day where # was the last digit in the Julian date. 
filel script 


edit] scnpt 251001 
edit script 251002 
edit 1] script 251003 
edit] script 251004 


edit] script 251086 


The “edit1l script” Unix script was the most important script. It uncompressed and 
compressed the raw data files, ran the stream editor command, executed the Fortran 


program, invoked convertTool and created the new edited data files. 


edit 1 script 


uncompress $1.sol - uncompresses raw data files 

sed -f sed1script $1.sol > dat.edt - runs stream editor and creates data file 
a.out - executes Fortran program 

cp out.ecf $1.ecf - copies Fortran program output to new file 
convertTool <out.ecf> $1.j20 - invokes convertTool program 

rm dat.edt - removes intermediary data file 

rm out.ecf - removes intermediary data file 

compress $1.sol - compresses raw data file for storage 
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The stream editor command eliminated countless hours of hand editing for the raw 
GPS data files by providing an automated way of removing unwanted file header lines and 
troublesome colons and slashes. The script “sediscript” is called by the SED command 


and contains the editing commands that were desired to be performed. 


sedlscript 

/GPS/d - deletes file header line which contains the character string GPS 
5/:/ /g - replaces all colons with spaces 

s!/! !g - replaces all slashes with spaces 


р. File Editing with Fortran Program 

The Fortran program “thingps” removed the desired state vector elements 
from the edited input file, converted time to MET in seconds, transformed GPS UEN 
velocity data into WGS 84 coordinates, created the new data files, wrote file header 
information required to execute convertTool and thinned the GPS data files by a factor of 


ten. An example of the output from “thingps” is shown in Figure 3.5. 


thingps.f 

program thingps 

implicit none - declare variables 
double precision x,y,z,velx,vely,velz,velup,veleast,velnorth 
double precision gps 

integer yr,doy,hr,min,sec,mo,day 

integer year,dayoyear,month,dayomon,hour,minute,second 
integer tmission,secs 

integer 1,count,11 


open(5,file-'dat.edt',status-'old") - open edited data file 
rewind(5) 

open(6,file='out.ecf' status='unknown') - open new data file 
rewind(6) 


year=1995 - data required for convertTool 
dayoyear=250 header and MET conversion 
month=09 

dayomon=07 

hour=15 

minute=9 

second=0 

tmission=dayoyear*24*3600+hour*3600+minute*60+second 
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% 


100 
200 


% 


% 


write(*,1) 3 - convert Tool header information 
format(11) 

write(*,2) year,month_dayomon,hour,minute,second 

format(14,5(5x,12)) 


count=0 
do 100 1=1,90000 
read (5,* ,end-200) gps,yr,doy,hr,min,sec x, y,z, - read desired data 
veleast,velnorth,velup 


secs=(doy*24*3600+hr*3600+min*60+sec)-tmission - MET conversion 
call uen(x,y,z,velup, veleast,velnorth, velx,vely, velz) 


i1 =1/10*10 - thin data files by factor of ten 
if (1.eq.11) then 

write(6,3) secs,x,y,z,velx, vely,velz 
end if 


count=countt | 
continue 
continue 
close(5) 
close(6) 
format(17,6(2x,f16.6)) 
end 


subroutine uen(x,y,z,velup, veleast, velnorth,velx,vely,velz) - velocity conversion 
implicit none 

double precision x,y,z,velup,veleast,velnorth 

double precision R,gamma,d,delta,velx,vely,velz 


Қас ту 722572005 - elements required for conversion 
gamma=atan2(y,x) (see UEN reference frame section) 
О (О ЫРА МЫСЫ 

delta=atan2(z,d) 


velx-cos(gamma)*cos(delta)*velup-sin(gamma)*veleast-cos(gamma) -velocity 
*sin(delta)*velnorth components 

vely=sin(gamma)*cos(delta)*velup+cos(gamma)*veleast- 
sin(gamma)*sin(delta)*velnorth 

velz=sin(delta)*velup+cos(delta)*velnorth 


return 
end 
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Figure 3.5. Output File from thingps.f. 
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& Reference Frame Transformation Using ConvertTool 

The convertTool program transformed the Fortran output file “out.ecf” 
from ECEF coordinates to J2000 coordinates producing a file with the extension “.j20”. 
The “out.ecf” header “3” corresponds to this conversion to J2000, and a “1” would 
correspond to a conversion from M50 coordinates to J2000 coordinates. ConvertTool 
also requires a start time, or epoch, in the file header when converting to or from ECEF 
coordinates. This start time was selected to be the time of launch in order to comply with 
the MET convention. Figure 3.6 shows an example output file from convertTool in J2000 


coordinates. 


3. Group GPS Files by Day 


After editing the raw GPS data files and creating the new files with time in MET 
and position and velocity in J2000 XYZ coordinates, the Unix script "daysscript" was 
utilized to group the data files by day. First, empty files were created for each day and 
labeled “gpsday#,” where # was equal to the last digit of the Julian date. Next, 
“daysscript” moved the empty files to a temporary file, “gps.j20,” for editing by the 
“append#script.” The “Append#scipt” invoked the script “dataappend” to append 
sequentially all the files for a flight day to create a single file. Finally, the edited files were 
moved back to the appropriate “gpsday#’ file. 


Daysscript append 1 script 
mv gpsday0.j20 gps.j20 dataappend 251001 
appendOscript dataappend 251002 
mv gps.j20 gpsday0.j20 dataappend 251003 
mv gpsday1.j20 gps.j20 | 
append 1script 
mv gps.j20 gpsday1.j20 | 
. dataappend 251086 
mv gpsday60.j20 gps.j20 dataappend 
append60scnpt ср $1.j20 data.j20 
mv gps.j20 gpsday60.j20 cat data.j20 >> gps.j20 
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Figure 3.6. ConvertTool Output File. 


30 


4. Edit Propagated Data File 


The propagated data file required much less manipulation. A Fortran program, 
“edtprop,” was used to add the header for convertTool, convert the input file times to 
MET in seconds and convert the data from English units to metric units. A copy of the 
program is in Appendix A. ConvertTool was then used to transform the file from M50 
coordinates to J2000 coordinates. Finally, the file of over 16,000 state vectors was hand 


edited to break it up by flight day. 


5. Create STK Vehicle Files 


STK Vehicle files for each flight day were created using the Shuttle STK vehicle 
library file as a guide. A vehicle file was created for each day from each data source. GPS 
ephemeris data and the propagated file ephemeris data were each inserted into a model 
Shuttle vehicle file in ECI J2000 coordinates. A scenario was created for every flight day 
which included a vehicle generated from GPS ephemeris data and one generated from the 
propagated ephemeris data. An example vehicle file is shown in Appendix C. The 
scenarios were animated so that both vehicles could be viewed simultaneously to better 


visualize the similarities and differences between the tracks. 


6. Matlab Comparison 


Matlab was employed for a detailed comparison of the GPS and propagated data 
files for each day. The Matlab M-file, “diffday#.m” was executed to plot J2000 XYZ 
position and velocity vector components and the position and velocity vector magnitudes 
for both data sets. A plot of the data points in three dimensions was also obtained using 
“orbplt#.m” to assist in visualizing the orbit that was generated by the GPS data. 
Appendices D and E contain examples of these M-files. 

Since the propagated data and GPS data could not be compared point for point, it 
was difficult to obtain hard numbers for comparison of GPS and reference track data. In 


an effort to match the data sets point for point, the process of editing the GPS data files 
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for days 251 and 252 was repeated with the exception of thinning the data by a factor of 
ten. The NAVG-11 file was hand edited to remove headers and unnecessary columns, 
and then run through the stream editor to remove colons. The Fortran program “edtsv” 
shown in Appendix B was then executed to convert state vector times to MET and 
English units to metric units. ConvertTool was invoked to transform the file from M50 
coordinates to J2000 coordinates. GPS and NAVG-11 state vectors were then hand 
matched for days 251 and 252. A Matlab M-file was used to compare position and 
velocity vector components, position and velocity vector magnitudes, and Root Mean 
Square (RMS) differences between the GPS data and the NAVG-11 data. A copy of the 
Matlab script “diffmat#.m” 1s shown is Appendix F. 
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IV. RESULTS 


GPS data was available for every flight day of STS-69; however, GPS data was 
not available for the entire duration of each flight day. The percentage of each day for 
which GPS data was available for comparison with the propagated data varied greatly 
from day to day. Day 258 was covered most sparsely by GPS with data being available 
only 12 pecent of the day. In contrast, GPS data was available for 93 percent of day 255. 
The average percentage of the day during the mission for which GPS data was available 


was 65 percent. These results are shown in Figure 4.1. 
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Figure 4.1. Percentage of each day for which GPS data was available. 
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A. MATLAB COMPARISON RESULTS 


1. Day 251 


Day 251 data was typical of the results achieved during the mission. Figure 4.2 at 
the end of the section displays how dense the GPS data is relative to the propagated file 


reference track. Of note is the fact that the GPS data, displayed as asterisks, has already 
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been thinned by a factor of ten. The propagated file is derived from the original state 
vectors which come from several sources as shown in Figure 4.3. The sources displayed 
trends when compared with the GPS data, and these trends appeared to be diurnal in some 


cases. The list of tracking source identifiers 1s shown in Table 4.1. 


Source Identifier 


East TDRS 
West TDRS 







Ascension Island 







Bermuda 





Kwajalein Island 


Wallops Island 





Table 4.1. NASA Shuttle Tracking Sources. 


A plot of the position vector XYZ components in J2000 coordinates for day 251 
shown in Figure 4.4 indicates that the GPS data represented by the heavy line follows very 
closely with the reference track represented by the thinner line. A two hour gap in GPS 
data is seen. This is typical of the gaps in GPS data that were observed every day. A plot 
of the velocity vector in XYZ components in J2000 coordinates in Figure 4.5 also shows 
this gap in GPS data and some spurious GPS data points. Both plots display sinusoidal 
patterns for each of the components which are associated with an inclined circular orbit. 

The position vector magnitude is equivalent to the spacecraft’s distance from the 
center of the Earth, or radius of its orbit, as shown in figure 4.6. Plotting of the position 
vector magnitude in Figure 4.7 shows that the radius is fairly constant as would be 
expected for the Shuttle's circular orbit. This plot of vehicle altitude was very effective in 
indicating when vehicle maneuvers occurred. The velocity vector magnitude plot also 
indicates that the velocity vector was quite constant as was expected for a circular orbit. 
Both the position and velocity vector magnitude plots show that spurious GPS data points 


were present. 
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The three dimensional plot of the GPS data for day 251 in J2000 coordinates 
shown in Figure 4.8 displays how well the GPS data captured the Shuttle’s circular orbit. 
The fact that the orbit is inclined, and that the GPS data is extremely dense is reflected in 
this plot. Erroneous data points are also visible. 

Point for point comparison of the GPS position vector magnitudes with the 
NAVG-11 position vector magnitudes in Figure 4.9 shows an altitude difference of +150 
m. Figure 4.10 shows that the velocity vector magnitude difference between the two 
data sets was generally + 1 m/s and did not exceed 4 m/s. The RMS position difference 
between the GPS data and NAVG-11 data displayed in Figure 4.11 was primarily between 
1400 and 1550 m, and the RMS velocity difference displayed in Figure 4.12 did not 
exceed 9 m/s. 


2: Typical Data 


The position and velocity vector magnitude plots for day 252 shown in Figure 
4.13 display a sinusoidal pattern associated with what appears to be an elliptical orbit. 
Closer examination of the position vector magnitude plot scale shows that the peak to 
peak difference associated with apogee and perigee is approximately 5000 meters. Gaps 
in GPS data and spurious data points are also experienced. Day 255 position and velocity 
vector magnitude plots displayed in Figure 4.14 are similar to those for day 252; however, 
the errant data points are more severe. 

GPS data for day 258 was very sparse with only 12 percent of the flight day 
covered by GPS data. The sparse amount of data allowed a detailed look at both data sets 
as shown in Figure 4.15 and displayed the imperfect nature of an actual vehicle orbiting 
the Earth. Although there were a few spurious data points, the GPS data followed the 
reference track so closely that it was difficult to discern the GPS data from the reference 


track. 
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de Maneuvers 


The position vector magnitude plot proved extremely useful in detecting Shuttle 
maneuvers. The day 253 position vector magnitude plot shown in Figure 4.16 displays a 
gap in GPS data just prior to the maneuver. A small gap in GPS data is experienced at the 
beginning of the maneuver which involves an altitude change of 10,000 m; however, the 
GPS track follows the reference track through the rest of the maneuver. Figure 4.17 also 
displays an example of GPS tracking through a maneuver. This maneuver is more severe 
involving a change in altitude of 60,000 m. A loss of GPS data is experienced in the 
middle of the maneuver with GPS data being reacquired before completion of the 


maneuver. 


4. Errant GPS Data 


There are several examples of errant GPS data over the course of the eleven flight 
days of STS-69. One of the best examples occurs on day 250 and is shown in Figure 4.18 
and Figure 4.19. Figure 4.18 is a three dimensional plot of the GPS data for day 250. 
This plot shows a segment of GPS data that fails to conform to the Shuttle's circular orbit. 
This period of bad data can be clearly identified in the position vector component plot 
shown in Figure 4.19. At the start of the day's data the GPS data does not follow the 
reference track. After a period of approximately 1300 seconds, the GPS data changes 
drastically and matches the reference track. 

GPS data for day 259 has an interesting anomaly that is exhibited in the position 
vector magnitude plot in Figure 4.20. A segment of GPS data appears to be disjointed 
from the rest of the data. In the three dimensional plot of Figure 4.21 this period of errant 
data appears as a separate orbit at a different inclination. Several spurious data points are 
also present, particularly in the velocity vector magnitude plot. 

The position and velocity vector magnitude plots presented the most dramatic 
results of errant data. The Day 260 position and velocity vector magnitude plot in Figure 
4.22 is an example of this observation. The last 5000 seconds of position and velocity 
data for the Day 260 GPS data varies significantly from the reference track. This variance 
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in the data includes as much as a 20 km difference in altitude and a 1.7 km/s difference in 


velocity. 
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Figure 4.2. Density of GPS Data. 
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Figure 4.3. NAVG-11 Data Sources. 
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Figure 4.4. Day 251 Position Vector Components in J2000 Coordinates. 
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Figure 4.5. Day 251 Velocity Vector Components in J2000 Coordinates. 
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Figure 4.6. Position Vector Magnitude. 
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Figure 4.7. Day 251 Position and Velocity Vector Magnitude. 
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Figure 4.8. Day 251 GPS Orbit in J2000 Coordinates. 
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Figure 4.9. Day 251 Position Vector Magnitude Difference. 
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Figure 4.10. Day 251 Velocity Vector Magnitude Difference. 
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Figure 4.11. Day 251 Position RMS Difference. 
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Figure 4.12. Day 251 Velocity RMS Difference. 
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Figure 4.13. Day 252 Position and Velocity Vector Magnitude. 
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Figure 4.14. Day 255 Position and Velocity Vector Magnitude. 
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Figure 4.15. Day 258 Position and Velocity Vector Magnitude. 
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Figure 4.16. Day 253 Position and Velocity Vector Magnitude. 
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Figure 4.17. Day 257 Position and Velocity Vector Magnitude. 
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Figure 4.18. Day 250 GPS Orbit in J2000 Coordinates. 
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Figure 4.19. Day 250 Position and Velocity Vector Magnitude. 
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Figure 4.20. Day 259 Position and Velocity Vector Magnitude. 
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Figure 4.21. Day 259 GPS Orbit in J2000 Coordinates. 
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Figure 4.22. Day 260 Position and Velocity Vector Magnitude. 
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B. STK RESULTS 


STK proved to be a valuable tool for visualizing the differences between the GPS 
vehicle tracks and the propagated vehicle tracks. The GPS vehicle 15 labeled 
“STS69gpsday#,” and the propagated file vehicle is labeled “STS69propday#” where # 
corresponds to the last digit of the Julian date. The circle surrounding the vehicle 
represents the vehicle's field of view as seen by a horizon sensor mounted on the Shuttle. 
This field of view circle was helpful in discerning changes in altitude because it expanded 
as the vehicle’s altitude increased. 

The ground tracks displayed in Figure 4.23 show an example of the extreme 
differences that were occasionally encountered between the two data sets. This example 
occurs on day 250 and provides another look at the discrepancy observed in the Matlab 
plots for this day shown in Figure 4.18 and Figure 4.19. The two tracks eventually did 
coincide as shown in Figure 4.24. When the vehicle tracks coincided in this manner it was 
difficult to distinguish between them because they were plotted directly on top of each 
Other. 

Figure 4.25 is more typical of the STK plots obtained. This plot shows the entire 
ground track for day 251. The GPS and propagated ground tracks overlay directly on top 
of each other, and the differences between the two vehicles are indistinguishable because 
they are also plotted on top of each other. Figure 4.26 displays a variance in the data for 
day 251 showing the divergent tracks of the two vehicles. Each day’s data analysis 
showed at least one anomaly similar to this one. The perspective shot in Figure 4.27 
provides a better view of the difference between the two tracks. Figure 4.28 shows a 
similar error for day 253. The expanded field of view for the GPS vehicle indicates that 
there is a difference in altitude between the two vehicles. 

A perspective view is shown in Figure 4.29 and displays how well the vehicles and 
orbits for day 255 match each other. Figure 4.30 is a close-up of the same picture and 
exhibits how closely the les and two sets of orbits overlay on top of each other. 
Figure 4.31 further zooms in and reveals that the vehicles are separated, with the GPS 


track trailing the propagated track by a small distance. 
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Several possibilities exist as to why the GPS data and reference track data do not 
match directly. One possibility is that the error is caused by the transformation between 
reference frames that does not take into account pole wander or the irregular rotation rate 
of the Earth. Another possibility is a spurious data point from the GPS receiver. While 
these factors may contribute to the difference between the data sets, the largest source of 
error is most probably the fact that the files for the GPS track and the files for the 
reference track each had times that were rounded to the nearest whole second. With the 
Shuttle traveling at over 7 km/s this can result in considerable position differences between 
the two data sets particularly in the downtrack direction. 

Figure 4.32 shows a GPS position data point at the center of a box whose sides 
were determined based on the velocity of the Shuttle in the X, Y and Z directions for a 
period of one second at a particular instant on day 251. This error box is approximately 
6000 m long in the X direction, 2500 m wide in the Y direction and 1200 m high in the Z 
direction. As expected, the reference track data point for the same instant fell within the 
box. Figure 4.33 shows a plot of the differences between the GPS data and the NAVG-11 
state vectors for day 251. The majority of the differences between the data sets lies within 
the box which encloses the error which may have been induced by a full second of 
timescale inaccuracy. 

The most extreme example of divergence between the GPS and propagated tracks 
occurred on day 256 and is shown in Figure 4.34. In this plot the perspective view is from 
150,000 km above the Earth and the small disc at the center is the Earth. Figure 4.35 is 
more indicative of the differences between vehicle tracks that were encountered. Typically 
an error in which the GPS track started to lead or lag the propagated vehicle would occur 
once per day. 

The ground track plot for day 258 shown in Figure 4.36 reflects the sparse data for 
this flight day. The perspective view in Figure 4.37 shows how closely the two vehicle 
tracks match. The close-up view in Figure 4.38 also shows that the GPS track coincides 
with the propagated track. Figure 4.39 reveals that the two tracks are nght on top of each 
other. The GPS orbit in this plot appears very smooth while the propagated orbit appears 
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to be segmented and not as smooth as the GPS orbit. This is expected since the GPS orbit 
is composed of considerably more data points. 

The plot for day 259 shown in Figure 4.40 indicates some discrepancies between 
the GPS and propagated vehicle ground tracks. It shows two independent vehicles which 
indicates that the vehicle tracks diverge. Figure 4.41 further shows that the GPS data has 
produced a vehicle in a different orbit at a different inclination. This plot matches the 


Matlab three dimensional GPS orbit plot for this day shown in Figure 4.21. 
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Figure 4.23. Day 250 STK Divergent Shuttle Vehicle Tracks. 
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Figure 4.24. Day 250 STK Coincident Shuttle Vehicle Tracks. 
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Figure 4.25. Day 251 STK Coincident Shuttle Vehicle Tracks. 
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Figure 4.26. Day 251 STK Ёс Shuttle Vehicle Tracks. 
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Figure 4.27. Day 251 STK Perspective View of Divergent Shuttle Vehicle Tracks. 
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Figure 4.28. Day 253 STK Divergent Shuttle Vehicle Tracks. 
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Figure 4.29. Day 255 STK Perspective View of Coincident Shuttle Vehicle Tracks. 
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Figure 4.30. Day 255 Perspective View Close-up of Shuttle Vehicle Tracks 
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Figure 4.31. Day 255 STK Perspective View Close-up of Shuttle Vehicle Tracks. 
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Day 251 GPS and Reference Pasitions in J2000 Coordinates 
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Figure 4.32. Day 251 GPS and Reference Positions in J2000 Coordinates. 
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Figure 4.33. Day 251 Error Box and Position Differences in J2000 Coordinates. 
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Figure 4.34. Day 256 STK Perspective View of Divergent Shuttle Vehicle Tracks. 
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Figure 4.35. Day 257 STK Perspective View of Divergent Shuttle Vehicle Tracks. 
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Figure 4.36. Day 258 STK Coincident Shuttle Vehicle Tracks. 
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Figure 4.37. Day 258 STK Perspective View of Coincident Shuttle Vehicle Tracks. 
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Figure 4.38. Day 258 STK Perspective View Close-up of Shuttle Vehicle Tracks. 
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Figure 4.39. Day 258 STK Perspective View Close-up of Shuttle Vehicle Tracks. 
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Figure 4.40. Day 259 STK Divergent Shuttle Vehicle Tracks. 
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Figure 4.41. Day 259 STK Perspective View of Divergent Shuttle Vehicle Tracks. 
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C. SHUTTLE UEN COMPARISON RESULTS 


1. UEN transformation 


In an effort to better classify the state vector errors encountered as errors in the 
crosstrack direction, in the downtrack direction, or in altitude, the GPS and NAVG-11 
state vectors with matching times were transformed from J2000 coordinates to the 
Shuttle’s UEN frame. In order to perform the transformation from the J2000 reference 
frame to the Shuttle UEN reference frame as shown in Figure 4.42, two Euler angle 


rotations must be performed. A rotation about the Z-axis by a value of y followed by a 


rotation about the Y-axis by -6 are required as shown in Equations 4.1, 4.2 and 4.3. 


Z 


CROSSTRACK 





Figure 4.42. Up-East-North Reference Frame. 


9 


U X 
Е| -С,-9С0)Ү (4.1) 
Z 


Shuttle J 2000 


а ô) 0 -sin(-ö) || cos(y) sın(y) 0 X 
1 0 —sin(y) cos(y) 0 Ё (4.2) 
мэн m ô) 0 со%-д9) 0 0 1 2122 
cos(d)cos(y) cos(d)sin(y)  sin(d) | | X 
E (4.3) 


U 
E = sin(6) cos(y) 0 
М - біп(д)сов(у) - віп(д)зіп(у) соя(д)| |2 |, 


Shuttle 


2: UEN results 


This transformation was implemented using Matlab and is available in Appendix G. 
The following results shown in Tables 4.2 and 4.3 and Figures 4.43, 4.44, 4.45, 4.49, 4.50 
and 4.51 at the end of the chapter were obtained from the day 251 and day 252 state 


vector position differences. 


їг oS Vertical (m) Downtrack (m) Crosstrack (m) 


Average 


Median 


Standard Deviation 
Table 4.2. Day 251 State Vector Position differences. 


ШИ = Vertical (m) Downtrack (m) Crosstrack (m) 


Average 





Median 


Standard Deviation 
Table 4.3. Day 252 State Vector Position differences. 
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The vertical position differences shown in Figure 4.43 and 4.49 were less than the 
expected one sigma value of 75 m for SPS. The values obtained for vertical position 
difference matched the values obtained from the position vector magnitude (altitude) 
difference plot when the data sets were analyzed in the J2000 reference frame. These 
values are within the expected SPS 156 m vertical positioning accuracy. 

Although the average downtrack position difference was fairly small for both days, 
the extremely large standard deviation indicates that significant downtrack errors were 
encountered. When traveling at 7 km/s a downtrack error of 1400 m corresponds to a 0.2 
s error in time. It was expected that a constant position error in the downtrack direction 
would be observed. The downtrack position differences displayed in Figure 4.44 and 
Figure 4.50 varied approximately +1400 m, and no constant offset was discovered. These 
results far exceeded the SPS one sigma horizontal positioning error of 50 m. 

Crosstrack errors were all less than 1x10% m for both days as shown in Figure 4.45 
and Figure 4.51. This was beyond the precision of the data, and makes the differences 
essentially zero. These results were unexpected and did not match the one sigma 
horizontal position error of 50 m for SPS. Due to unidentified sources of error, no firm 


conclusions can be drawn; however, the downtrack and crosstrack errors seem to indicate 


that most of the position error is in the downtrack direction. 


БЕ Vertical (m/s) Downtrack (m/s) Crosstrack (m/s) 


Average 


Median 


Standard Deviation 





Table 4.4. Day 251 State Vector Velocity differences. 


EU — же Vertical (m/s) Downtrack (m/s) Crosstrack (m/s) 
1.13 | 


Average 


Median 


Standard Deviation 





Table 4.5. Day 252 State Vector Velocity differences. 
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The velocity results from the comparison of state vectors for day 251 and day 252 


are shown in Tables 4.4 and 4.5 and Figures 4.46, 4.47, 4.48, 4.52, 4.53 and 4.54. Both 


days’ results are very similar. In general the average values for vertical and downtrack 


velocity differences exceeded the expected SPS one sigma velocity accuracy of 0.5 m/s, 


and the median values were close to this value. The crosstrack velocity data is well within 


the expected SPS one sigma velocity accuracy. 
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Figure 4.43. Day 251 Vertical Position Difference. 
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Downtrack Position Difference (m) 
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Figure 4.44. Day 251 Downtrack Position Difference. 
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Crosstrack Position Difference (m) 
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Figure 4.45. Day 251 Crosstrack Position Difference. 
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Vetical Velocity Difference (m/s) 
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Figure 4.46. Day 251 Vertical Velocity Difference. 
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Figure 4.47. Day 251 Downtrack Velocity Difference. 
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Figure 4.48. Day 251 Crosstrack Velocity Difference. 
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Vetical Position Difference (m) 
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Figure 4.49. Day 252 Vertical Position Difference. 
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Figure 4.50. Day 252 Downtrack Position Difference. 
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Figure 4.51. Day 252 Crosstrack Position Difference. 
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Vetical Velocity Difference (m/s) 
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Figure 4.52. Day 252 Vertical Velocity Difference. 
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Figure 4.53. Day 252 Downtrack Velocity Difference. 
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Figure 4.54. Day 252 Crosstrack Velocity Difference. 
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Ve CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


GPS navigation solutions were available for approximately 65 percent of the STS- 
69 mission, and they generally coincided with the reference track. Based on vehicle track 
visualization using STK and based on conversion from J2000 coordinates to a spacecraft 
local UEN reference frame, differences between the GPS data and the reference data 
appear to have occurred predominantly in the downtrack direction. The fact that the two 
data sets were rounded to the nearest whole second probably impacted the downtrack 
differences between them significantly. 

State vector differences between the GPS navigation solutions obtained using the 
Standard Positioning Service and the reference NAVG-11 state vectors for day 251 and 
252 produced RMS position differences between the data sets of about 1500 m. One 
sigma position accuracy of 54 m in the vertical direction and approximately 1400 m in the 
downtrack direction were experienced. Velocity vector magnitude differences during this 
period were generally + 1m/s, with a RMS velocity difference of less than 9 m/s. One 
sigma velocity accuracies of approximately 4.2 m/s in the vertical direction, 2.3 m/s in the 
downtrack direction and 1.5 m/s in the crosstrack direction were experienced. A firm 
conclusion regarding GPS accuracies could not be drawn because all sources of error 
were not identified. 

GPS position and velocity data was generally very good; however, spurious data 
points were present. The errors included in these data points were sometimes very 
significant, and outages in GPS data of approximately two hours appeared to occur at 
least once per day. Despite these errors, GPS appeared to be effective in producing good 
state vector data even during vehicle maneuvers. 

Based on these results GPS appears to be an excellent navigation source for 
providing Shuttle state vector information. Considering that this data can be obtained 
real-time and still match the reference so closely, real-time GPS state vector data appears 


even more valuable. An important caveat must be applied to this very valuable GPS data: 


25) 


another navigation source such as INS must be present to provide a check against 


spurious data points and periods of outage. 


B. RECOMMENDATIONS 


Analysis of state vector data from a Precise Positioning Service GPS receiver 
should be conducted. PPS state vector data should provide considerably better position 
and velocity data which is required during rendezvous and proximity operations. STS-79 
will provide an opportunity to assess PPS GPS receiver navigation solutions. 

Before integration with the Shuttle navigation suite, experimentation with a simple 
filter for the GPS receiver navigation solutions is required. Erroneous GPS state vector 
data was sometimes quite severe. A process such as a Kalman filter is required to remove 
navigation solutions outside a specified tolerance. 

A standard library of code in a commonly accepted language such as Fortran 
should be adopted for frequently required astronautical functions such as conversion 
between reference frames. In particular, a simple uniform routine for transformation 
between the J2000 ECI frame and the WGS-84 ECEF frame using general precession 
according to IAU-1976 and general nutation according to IAU-1980 without considering 
pole wander or the irregular rotation of the Earth is required. Algorithms which allow real 
time processing of data such as GPS navigation solutions should be emphasized. In 
contrast, current references emphasize the use of techniques which require post-flight 
processing due to their requirement of celestial observations. 

GPS system time should be adopted by all tracking sources. Ambiguity caused by 
UTC being rounded to the nearest whole second by the GPS receiver and tracking sources 
led to significant errors, particularly in the downtrack direction. Since GPS system time 
can be obtained to a very high accuracy by both the Shuttle receiver and a receiver at the 
tracking site, GPS time would provide a valuable timescale addition to all tracking data 


sets. 
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APPENDIX A. FORTRAN PROGRAM FOR EDITING PROPAGATED FILE 
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wr sous IYI 


100 





100 
200 


program edtprop 


implicit none - declare variables 
double precision hours, xft, yft,zft, vxft, vyft, vzft 

double precision secs,x,y,Z,VX,VY,VZ 

integer 1,count 


open(5,file='prop.out',status='old') - open data file 

rewind(5) 

open(6,file-'prop.m50',status-'unknown') - open new data file 

rewind(6) 

write(6,*) 1 - convertTool header to 
transform to J2000 

count=0 


do 100 i=1,90000 
read (5,* end=200) hours,xft, yft,zft_vxft_vyft vzft -read data 


secs=hours*3600.0 - covert MET hours to seconds 
x=xft*0.3048 - Convert feet to meters 
y=yft*0.3048 

z=zft*0.3048 

vx=vxft*0.3048 


vy=vyft*0.3048 
vz=vzit* 0.3048 


write(6,1) secs,x,y,Z,VX,VY,VZ - write to new data file 
count=count+1 
continue 


continue 


close(5) 
close(6) 


format(f16.6,6(2x,f16.6)) 
end 
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APPENDIX B. FORTRAN PROGRAM FOR EDITING NAVG-11 STATE 
VECTOR FILE 
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9/0 


program edtsv 


implicit none - declare variables 
double precision station, xft, yft,zft, vxft,vyft, vzft 

double precision x,y,z,vx,vy,vz 

integer sv,doy,hr,min,sec 

integer year,dayoyear,month,dayomon,hour,minute,second 
integer tmission,secs 

integer i,count 


open(5,file-'sv.edt' status-'old') - open data file 
rewind(5) 

open(6,file='sv.m50',status='unknown') - open new data file 
rewind(6) 


year=1995 - data required for MET 
dayoyear=250 conversion 

month=09 

dayomon=07 

hour=15 

minute=9 

second=0 


tmission=dayoyear*24*3600+hour*3600+minute*60+second - MET start 


write(*,1) 1 - convertTool header to 
format(11) transform to J2000 
count=0 


do 100 i=1,90000 
read (5,*,end=200) station,sv,doy,hr,min,sec, - read desired data 
xft,yft,zft,vxft,vyft,vzft 


secs=(doy*24*3600+hr*3600+min*60+sec)-tmission - MET conversion 


x=xft*0.3048 - convert feet to meters 


y=yft*0.3048 
2-28%0.3048 


vx=vxft*0.3048 
vy=vyft*0.3048 


vz=vzft* 0.3048 


write(6,2) secs,x,y,z, vx,vy,vz - write to new data file 


count=count+1 
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100 
200 


continue 
continue 


close(5) 
close(6) 


format(i7,6(2x,f16.6)) 
end 
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APPENDIX C. STK VEHICLE FILE 
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stk.v.2.0 


BEGIN Vehicle 


Name 


BEGIN VehiclePath 


BEGIN J4Perturbation 


EphemEpoch 


StartTime 
StopTime 
SemiMajorAxis 
Altitude OfApogee 
AltitudeOfPerigee 
RadiusOfApogee 
RadiusOfPerigee 
PeriodOfOrbit 
MeanMotion 
Inclination 
Eccentricity 
ArgOfPerigee 
LongAscenNode 
Right Ascension 
TrueAnomaly 
MeanAnomaly 


TimePast AscenNode 


TimePastPerigee 
TimeStep 

OrbitGranularity 
NumberOfPasses 


SizeShapeT ype 
AscenNodeType 
AnomalyType 
EllipseType 


END J4Perturbation 


BEGIN PassDefn 


Break 
Latitude 
DisplayFlag 


STS69gpsday 1 


7 Sep 1995 15:09:00.000000 


7 Sep 1995 15:09:00.00 
19 Sep 1995 00:00:00.00 
6748537.000000 
370400.000000 
370400.000000 
6748537.000000 
6748537.000000 
5517.284751 
0.001139 
28.500000 
0.000000000000 
0.000000 
-151.000000 
249.528639 
0.000000 
0.000000 
0.000000 
0.000000 
60.000000 
0.000000 
0 


АррРепАН 
Longitude 
TrueAnomaly 
Osculating 


Ascending 
0.000000 
Both 
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FirstPass 1 


RangeFirstPass 0 
RangeLastPass 2147483647 
DisplayScheme AllEphemeris 
ScenarioEpoch ] Nov 1992 00:00:00.0 
Passes 
END Passes 

END PassDefn 

END VehiclePath 

BEGIN Ephemeris 

NumberOfEphemerisPoints 17342 

ScenarioEpoch 1 Хоу 1992 00:00:00.0 

EphemerisEcfTimePosVel 


89944473 .000000 -629041 1.999970 -10669.639545 2444894 749963 -944 961651 - 
7710.190409 -2382.491604 

89944485 000000 -6300850 499989 -92193.726556 2416095.999986 -847.046094 - 
7710.033323 -2420.407165 

89944522.000000 -6325375.499997 -33 1401.656250 2329124.749954 -591.753766 - 
7704.550020 -2514.460050 


90020994.000000 -5868722.494942 883252.418377 3213828.661367 -1541.348168 - 
7851.489098 -89.026636 
90020994.000000 -5868722.494942 883252.418377 3213828.661367 -1541.348168 - 
7851.489098 -89.026636 
90020994.000000 -5868722.494942 883252.418377 3213828.661367 -1541.348168 - 
7851.489098 -89.026636 
END Ephemeris 
BEGIN Attitude 

BEGIN ECFVVLH 

AZIMUTH 0.000000 
END ECFVVLH 


END Attitude 
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BEGIN Swath 


ElevationAngle 
HalfAngle 
RepType 


END Swath 
BEGIN Constraints 


ConstraintMask 
MinRange 
MaxRange 
MinAzimuthAngle 
MaxAzimuthAngle 
MinElevationAngle 


MaxElevationAngle 


MinGrazingAngle 
MaxGrazingAngle 


MinGrazing Altitude 
MaxGrazing Altitude 


0.000000 
0.000000 
NoSwath 


0 

0.000000 
1000000000.000000 
-90.000000 
180.000000 
-45.000000 
45.000000 
0.000000 
45.000000 
185200.000000 
1852000.000000 


MinGndElevationAngle 0.000000 
MaxGndElevationAngle 90.000000 
MinSunElevationAngle -90.000000 
MaxSunElevationAngle 90.000000 


LightingCondition 
END Constraints 
BEGIN Attributes 
MarkerColor 
GroundTrackColor 
SwathColor 
LineStyle 
MarkerStyle 
END Attributes 
BEGIN Graphics 

Inherit 


ShowLabel 
ShowGndTrack 


5 


BWANA On 


On 
On 
On 
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END Graphics 
BEGIN Extensions 


BEGIN Desc 
END Desc 


BEGIN Group 
END Group 


BEGIN Eclipse 
Penumbra O 
PenumbraColor 7 
0 
1 


PenumbraLineStyle 
PenumbraLineWidth 


Sunlight O 
SunlightColor 7 
SunlightLineStyle 0 
SunlightLineWidth 1 
Umbra Off 
UmbraColor 7 
UmbraLineStyle 0 
UmbraLineWidth 1 
END Eclipse 


BEGIN Exclusion 
END Exclusion 


END Extensions 

BEGIN SubObjects 

Class Sensor 
Horizon 

END Class 

END SubObjects 


END Vehicle 
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APPENDIX D. MATLAB M-FILE FOR COMPARING GPS AND 
PROPAGATED DATA 


113 
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о 


% Program: diffdayl.m 
% Jım Jones 
% 25 May 96 


e 


% This program displays x,y and z components of GPS and reference 
% track position and velocity data for one day of data. It also 

Vo displays the magnitude of the position vector, or altitude, 

Vo from the center of the earth and the magnitude of the velocity 

% vector. 


% Load the GPS data file for one day 


load gpsday1.j20 -ascii 


gpsin-gpsday1; 


% Create time, position component and velocity component vectors. 


gpssec=gpsin(:,1); 
gpsx=gpsin(:,2); 
gpsy=gpsin(:,3); 
gpsz=gpsin(:,4); 
gpsxdot-gpsin(:,5); 
gpsydot-gpsin(:,6); 
gpszdot-gpsin(:,7); 


% Load the reference track data file for one day 
load propday1.j20 -ascil 


propin=propdayl; 


% Create time, position component and velocity component vectors. 


propsec-propin(:,1); 
propx-propin(:,2); 
propy=propin(:,3); 
propz-propin(:,4); 
propxdot-propin(:,5); 
propydot=propin(:,6); 
propzdot-propin(:,7); 
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% Plot position components vs. time 


figure(1) 

subplot(3,1,1) 

plot(gpssec,gpsx,'w*' propsec,propx,'w') 
xlabel('Day 251 Mission Seconds") 
ylabel( X Position (m)') 

grid 


subplot(3,1,2) 
plot(gpssec,gpsy,'w*',propsec,propy,'w') 
xlabel('Day 251 Mission Seconds") 
ylabel('Y Position (m)') 

grid 


subplot(3,1,3) 

plot(gpssec,gpsz,'w*' propsec,propz,'w') 
xlabel('Day 251 Mission Seconds") 
ylabel(Z Position (m)') 

tol=-1; 

legend('GPS','Reference’,tol) 

grid 

orient tall 

Yoprint_onig -dgif8 251 pos. gif 

print 


% Plot velocity components vs. time 


figure(2) 

subplot(3,1,1) 

plot(gpssec,gpsxdot,'w*' propsec,propxdot,'w') 
xlabel(Day 251 Mission Seconds’) 

ylabel('X Velocity (m/s)) 

grid 

subplot(3,1,2) 

plot(gpssec,gpsydot,'w*' propsec,propydot,'w') 
xlabel('Day 251 Mission Seconds") 

ylabel("Y Velocity (m/s)') 

grid 
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subplot(3,1,3) 

plot(gpssec,gpszdot,'w*' propsec,propzdot,'w') 
xlabel('Day 251 Mission Seconds') 

ylabel(Z Velocity (m/s)') 

tol=-1; 

legend('GPS',Reference' tol) 

grid 

orient tall 

Yoprint_ong -dgif8 251 vel. gif 

print 


?o Calculate position vector magnitudes (altitude) 


gpsr-sqrt(gpsx.^2 + gpsy.”2 + gpsz.”2); 
propr-sqrt(propx.^2 ^ propy.^2 ^ propz.^2); 


?o Calculate velocity vector magnitude 


gpsv-sqrt(gpsxdot.^2 ^ gpsydot.^2 * gpszdot.^2); 
propv-sqrt(propxdot.^2 + propydot.*2 + propzdot.’2); 


% Plot position and velocity vector magnitudes 


figure(3) 

subplot(2, 1,1) 

plot(gpssec,gpsr,'w*' propsec,propr,'w') 

xlabel(‘Day 251 Mission Seconds’) 

ylabel(‘Position Vector Magnitude (Vehicle Altitude) (ту) 
tol=-1 

legend('GPS' Reference',tol) 

grid 


subplot(2,1,2) 

plot(gpssec,gpsv,'w*' propsec,propv,'w”) 
xlabel('Day 251 Mission Seconds) 
ylabel("Velocity Vector Magnitude (m/s)’) 
tol=-1 

legend('GPS', Reference',tol) 

grid 

orient tall 

%print_ong -dgif8 251mag.gif 

print 


% Plot close up of position data 
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figure(4) 

plot(gpssec,gpsx,'w*' propsec,propx,'wo') 
xlabel('Day 251 Mission Seconds") 

ylabel( X Position (m)') 

tol=-] 

legend('GPS','Reference’,tol) 

grid 

ах18(((5е4-20) (5е4-220) -0.57е7 -0.47е71) 
Yoprint_ong -dgif8 251 close. gif 

print 
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APPENDIX E. MATLAB M-FILE FOR PLOTTING GPS ORBIT IN 3-D 
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% Program: orbplti.m 
% Jim Jones 
% 25 May 96 


% This program creates a 3D plot in J2000 coordinates of GPS navigation 
o solutions for one day of flight data. 


= 


% Load GPS ephemeris file in X,Y,Z J2000 coordinates 
load gpsday0.}20 -ascii 

out = gpsday0; 

% Plot ephemeris from data file 


х1 = [0 1е7]; 
х2 = [0 0]; 
x3 = [00]; 


yl = [00]; 
у2 = [0 1е7]; 
уЗ = [0 0]; 


21 = [00]; 
22 = [0 0]; 
23 = [0 1е7]; 


figure(1) 
plot3(x1,x2,x3,'w-',y1,y2,y3,'w-',Z1,Z2,z3,'w-',out(:,2),out(:,3 ),out(:,4),'w.') 
title'GPS Orbit for Day 250 in J2000 Coordinates") 

xlabel("X Axis (т)) 

ylabel("Y Axis (m)') 

zlabel('Z Axis (m)') 

%print_orig -dgif8 2500rb.gif 

print 
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=| xd 


APPENDIX F. MATLAB M-FILE FOR COMPARING GPS AND NAVG-11 
STATE VECTORS WITH MATCHING TIMES 
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HDV AAA PAD IMA NICO НОР АН НИ 
AA у ИУ 
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% Program: diffmatl.m 
% Jim Jones 
% 25 May 96 


% This program plots differences between GPS navigation solution 
% data and Rockwell state vector data which share the same time. 
% Both data sets are in J2000 coordinates. 


% Load GPS data file 

load gpsmatch1.j20 -ascıı 

gpsin-gpsmatchl; 

% Create time, position component and velocity component vectors. 


gpssec=gpsin(:,1); 
gpsx-gpsin(:,2); 
gpsy-gpsin(:,3); 
gpsz=gpsin(:,4); 
gpsxdot=gpsin(:,5); 
gpsydot-gpsin(:,6); 
gpszdot-gpsin(:,7); 


% Load Rockwell state vector data 

load svmatch1.j20 -ascıı 

propin=svmatchl; 

% Create time, position component and velocity component vectors. 


propsec=propin(:,1); 
propx=propin(:,2); 
propy=propin(:,3); 
propz=propin(:,4); 
propxdot-propin(:,5); 
propydot=propin(:,6); 
propzdot-propin(:,7); 
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% Identify sources of Rockwell state vectors 
% (Sources corresponding to indices of matrix propin) 


9/0 Legend 
estr=[1 3 5 79 11 16 21 24 27 30 32 34 36]; Фо * 
wstr=[2 4 6 8 10 13 17 18 22 28 85 37]; % о 
bdqc=[12 15 20 26 31]; % x 
wirc=[14 19 25]; % + 
kmtc=[23 29 33]; % . 


% Plot position components vs. time 


figure(1) 

subplot(3,1,1) 

plot(gpssec,gpsx,'w*', propsec,propx, wo’) 
xlabel(‘Day 251 Mission Seconds’) 
ylabel('X Position (m)') 

grid 

subplot(3,1,2) 

plot(gpssec,gpsy,'w*' propsec,propy,'wo') 
xlabel(‘Day 251 Mission Seconds’) 
ylabel('Y Position (m)') 

grid 


subplot(3, 1,3) 
plot(gpssec,gpsz,'w*',propsec,propz, wo’) 
xlabel('Day 251 Mission Seconds”) 
ylabel('Z Position (m)') 

tol=-1; 

legend('GPS','Reference’,tol) 

grid 

orient tall 

Yprint_ong -dgif8 251posm. gif 

print 


% Plot velocity components vs. time 


figure(2) 

subplot(3,1,1) 

plot(gpssec,gpsxdot,'w*' propsec,propxdot,'wo') 
xlabel('Day 251 Mission Seconds”) 

ylabel(X Velocity (m/s)') 


grid 
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subplot(3,1,2) 

plot(gpssec,gpsydot,'w*' propsec,propydot,'wo') 
xlabel('Day 251 Mission Seconds) 

ylabel(Y Velocity (m/s)) 

опа 


subplot(3,1,3) 

plot(gpssec,gpszdot,'w*' propsec,propzdot,'wo') 
xlabel(‘Day 251 Mission Seconds’) 

ylabel('Z Velocity (m/s)') 

tol=-1; 

legend('GPS','Reference’,tol) 

end 

orient tall 

%print_orig -dgif8 251velm.gif 

print 


% Calculate position vector magnitudes (altitude) 


gpsr-sqrt(gpsx.^2 ^ gpsy.^2 + gpsz.^2); 
propr-sqrt(propx.^2 ^ propy.^2 + ргорг.72); 


% Calculate velocity vector magnitude 


gpsv-sqrt(gpsxdot.^2 + gpsydot.*2 ^ gpszdot.^2); 
propv-sqrt(propxdot.^2 + propydot.”2 + propzdot.”2); 


% Plot position vector and velocity magnitudes 


figure(3 ) 

subplot(2, 1,1) 

plot(gpssec,gpsr,'w*', propsec, propr, wo’) 

xlabel(‘Day 251 Mission Seconds’) 

ylabel(‘Position Vector Magnitude (Vehicle Altitude) (m)') 
grid 


subplot(2,1,2) 

plot(gpssec,gpsv,'w*' propsec,propv,'wo') 
xlabel('Day 251 Mission Seconds") 
ylabel('Velocity Vector Magnitude (m/s)’) 
tol=-1 

legend(GPS' Reference' tol) 

grid 

orient tall 

Vepnrnt orig -dgif8 251magm. gif 
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print 

% Calculate position component differences 
posdiffx=gpsx-propx; 

posdifty=gpsy-propy, 

posdiffz-gpsz-propz; 


% Plot position component differences 


% Data source legend 


% estr - * 

% wstr- O 

Фо Бадс - х 

% wlrc - + 

% kmtc- . 
figure(4) 


plot(propsec(estr),posdiffx(estr),'w*',... 
propsec(wstr),posdiffx(wstr),'wo',‚propsec(bdqc),posdiffx(bdac),'wx',... 
propsec(wlrc),posdiffx(wirc),'w+',propsec(kmtc),posdiffx(kmtc),'w.’) 

xlabel('Day 251 Mission Seconds 

ylabel('X Position Difference (m)') 

grid 

tol=-1; 

legend ESTR'," WSTR','BDQC', WLRC','KMTC'tol) 

Voprint orig -dgif8 251difxm.gif 

print 


figure(5) 

plot(propsec(estr),posdiffy(estr),'w*',... 
propsec(wstr),posdiffy(wstr),'wo',propsec(bdqc),posdiffy(bdqc),'wx",... 
propsec(wIrc),posdiffy(wirc),'w+',‚propsec(kmtc),posdiffy(kmtc),'w.') 

xlabel('Day 251 Mission Seconds?) 

ylabel('Y Position Difference (m)') 

grid 

tol=-1; 

legend(ESTR''WSTR''BDOQC''WLRC''KMTC', tol) 

%print_orlg -dgif8 251difym. gif 

print | 
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figure(6) 

plot(propsec(estr),posdiffz(estr),'w*',... 
propsec(wstr),posdiffz(wstr),'wo',propsec(bdqc), posdiffz(bdqc),'wx’,.... 
propsec(wlrc), posdiffz(wirc),'wt',propsec(kmtc), posdiffz(kmtc),'w.') 

xlabel("Day 251 Mission Seconds’) 

ylabel('Z Position Difference (т)) 

grid 

tol=-1; 

legend('ESTR' 'WSTR''BDQC''WLRC''KMTC' tol) 

%ргіпі orig -dgif8 251difzm.gif 

print 


% Calculate velocity component differences 


veldiffx=gpsxdot-propxdot; 


veldiffy-gpsydot-propydot; 
veldiffz-gpszdot-propzdot; 


% Plot velocity component differences 


figure(7) 

plot(propsec(estr),veldiffx(estr),'w*',... 
propsec(wstr),veldiffx(wstr),'wo',propsec(bdqc),veldiffx(bdgqc),'wx',... 
propsec(wirc),veldiffx(wlrc),'w+',propsec(kmtc),veldiffx(kmtc),'w.’) 

xlabel('Day 251 Mission Seconds) 

ylabel(X Velocity Difference (m/s)) 

grid 

tol=-1; 

legend('ESTR''WSTR''BDQC''WLRC''KMTC' tol) 

%print_orig -dgif8 251dfxdm.gif 

print 


figure(8) 

plot(propsec(estr),veldiffy(estr),'w*',... 
propsec(wstr),veldiffy(wstr),'wo',propsec(bdqc), veldiffy(bdqc),'wx’,... 
propsec(wirc),veldiffy(wirc),'w+',propsec(kmtc),veldiffy(kmtc),'w.') 

xlabel('Day 251 Mission Seconds) 

ylabel("Y Velocity Difference (m/s)') 

grid 

tol=-1; 

legend('ESTR''WSTR''BDOC''WLRC''KMTC', tol) 

%print_orig -dgif8 251dfydm.gif 

print 
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figure(9) 

plot(propsec(estr),veldiffz(estr),'w*',... 
propsec(wstr),veldiffz(wstr),'wo',propsec(bdqc), veldiffz(bdqc),'wx’.... 
propsec(wIrc),veldiffz(wIrc),'w+',propsec(kmtc),veldiffz(kmtc),'w.') 

xlabel('Day 251 Mission Seconds") 

ylabel(Z Velocity Difference (m/s)') 

grid 

tol=-1; 

legend(ESTR''WSTR''BDQC''WLRC''KMTC' tol) 

%print_ong -dg1f8 251dfzdm. gif 

print 


% Calculate position and velocity magnitude differences 
diffr-gpsr-propr; 

diffv-gpsv-propv; 

% Plot position and velocity magnitude differences 


figure(10) 

plot(propsec(estr),diffr(estr),'w*',... 
propsec(wstr),diffr(wstr),'wo',propsec(bdqc),diffr(bdqc),'wx',... 
propsec(wirc),diffr(wirc),'w+',propsec(kmtc),diffr(kmtc),'w.') 

xlabel(‘Day 251 Mission Seconds’) 

ylabel(‘Position Vector Magnitude (Altitude) Difference (m)') 

grid 

tol=-1; 

legend(ESTR' 'WSTR''BDQC''WLRC''KMTC' tol) 

Yoprint_ong -dgif8 251difrm. gif 

print 


figure(11) 

plot(propsec(estr),diffv(estr),'w*',... 
propsec(wstr),diffv(wstr),'wo',propsec(bdqc),diffv(bdqc),'wx',... 
propsec(wlrc),diffv(wlrc),'w+', propsec(kmtc),diffv(kmtc),'w.') 

xlabel(‘Day 251 Mission Seconds’) 

ylabel('Velocity Vector Magnitude Difference (m/s)') 

grid 

tol=-1; | 

legend(ESTR''WSTR''BDQC''WLRC''KMTC' tol) 

Veprint orig -dgif8 25 1difvm. gif 

print 
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% Calculate position and velocity RMS differences 
rmsr-sqrt(posdiffx.^2 -- posdiffty.^2 -- posdiffz.^2); 
rmsv-sqrt(veldiffx.^2 4- veldiffy.^2 -- veldiffy.^2); 
% Plot position and velocity RMS differences 


figure(12) 

plot(propsec(estr),rmsr(estr),'w*',... 
propsec(wstr),rmsr(wstr),'wo',‚propsec(bdge),rmsr(bdqc),'wx‘',... 
propsec(wirc),rmsr(wirc),'w+',propsec(kmtc),rmsr(kmtc),'w.’) 

xlabel('Day 251 Mission Seconds") 

ylabel('Position RMS Difference (m)') 

grid 

tol=-1; 

legend(ESTR''WSTR''BDOQC''WLRC''KMTC' tol) 

%print_orig -dgif8 251drmsm. gif 

print 


figure(13) 

plot(propsec(estr),rmsv(estr),'w*",... 
propsec(wstr),rmsv(wstr),'wo',propsec(bdqc),rmsv(bdqc),'wx’,... 
propsec(wlrc), rmsv(wlrc),'w+',propsec(kmtc),rmsv(kmtc),'w.’) 

xlabel('Day 251 Mission Seconds) 

ylabel( Velocity RMS Difference (m/s)') 

grid 

tol=-1; 

legendESTR''WSTR''BDQC''WLRC''KMTC' tol) 

%print_orig -dgif8 251dvmsm.gif 

print 
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APPENDIX G. MATLAB M-FILE FOR TRANSFORMING MATCHING GPS 
AND NAVG-11 STATE VECTORS TO UEN FRAME FOR COMPARISON 
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% Program: uen251.m 
% Jim Jones 
% 25 June 96 


% This program plots differences between GPS navigation solution 
% data and NAVG-11 state vector data which share the same time. 
9o Both data sets are in the Shuttle UEN frame. 


% Load the GPS data file for one day 

load gpsmatch1.j20 -ascii 

gpsin=gpsmatchl; 

% Create time, position component and velocity component vectors. 


gpssec-gpsin(:,1); 
gpsx=gpsin(:,2); 
epsy=gpsin(:,3); 
gpsz=gpsin(:,4); 
gpsxdot-gpsin(:,5); 
gpsydot-gpsin(:,6); 
gpszdot=gpsin(:,7); 


% Load the reference track data file for one day 

load svmatch1.j20 -ascıı 

svin=svmatchl; 

% Create time, position component and velocity component vectors. 
refsec=svin(:,1); 

refx=svin(:,2); 

refy=svin(:,3); 

refz=svin(: ‚4); 

refxdot=svin(:,5); 


refydot=svin(:,6); 
refzdot=svin(:,7); 
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% Identify sources of NAVG-11 state vectors 
% (Sources corresponding to indices of matrix svin) 


% Legend 
estr=[1 35 791116 21 24 27 30 32 34 36];% * 
wstr=[2 4681013 17 18 22 28]; % O 
bdqc=[12 15 20 26 31); % X 
wlrc-[14 19 25]; Por 
kmtc=[23 29 33]; v e 


for 1=1:length(gpsx) 


% 


% 


% 


% 


Elements required to calculate DCM for GPS data 


gpsR-sqrt(gpsx(1)^2*gpsy(1)^2*gpsz(1)^2); 
gpsgamma-atan2(gpsy(1),gpsx(1)); 
gpsd=sqrt(gpsx(1)*2+gpsy(1)"2); 
gpsdelta=atan2(gpsz(i),gpsd); 


DCM components for GPS data 

gpsdcm1 1=cos(gpsdelta)*cos(gpsgamma); 
gpsdcm12-cos(gpsdelta)*sin(gpsgamma); 
gpsdcm13-sin(gpsdelta); 
gpsdcm21=sin(gpsdelta); 
gpsdcm22=cos(gpsgamma); 

gpsdcm23=0; 

gpsdcm3 1=-sin(gpsdelta)*cos(gpsgamma); 
gpsdcm32--sin(gpsdelta)*sin(gpsgamma); 
epsdcm33=cos(gpsdelta); 

DCM rows for GPS data 
gpsdcm1=[gpsdcm11 gpsdcm12 gpsdcm13 ]; 
gpsdcm2-[gpsdcm21 gpsdcm22 gpsdcm23]; 
gpsdcm3=[gpsdcm31 gpsdcm32 gpsdcm33]; 
3x3 DCM for GPS data 


gpsdcm=[gpsdcm1; gpsdcm?2; gpsdcm3]; 
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% 


% 


% 


% 


% 


% 


% 


Create GPS position and velocity vectors 


gpspos-[gpsx(1); gpsy(i); gpsz(i)]; 
gpsvel-[gpsxdot(1); gpsydot(1); gpszdot(1)]; 


Perform transformation from J2000 coordinates to UEN 


gpsuenp=gpsdcm* gpspos; 
gpsuenv-gpsdcm* gpsvel; 


Elements required to calculate DCM for reference data 


refR-sqrt(refx(1)^2-refy(1)^2-refz(1)^2); 
refgamma=atan2(refy(1),refx(1)); 
refd-sqrt(refx(1)^2-*refy(1)^2); 
refdelta=atan2(refz(i),refd); 


DCM components for reference data 
refdcm1 1=cos(refdelta)* cos(refgamma); 
refdcm12-cos(refdelta)*sin(refgamma); 
refdcm13-sin(refdelta); 

refdcm2 1=sin(refdelta); 
refdcm22-cos(refgamma); 

refdcm23-0; 

refdcm3 1=-sin(refdelta)*cos(refgamma); 
refdcm32=-sin(refdelta)*sin(refgamma); 
refdcm33=cos(refdelta); 

DCM rows for reference data 
refdcm1=[refdcm11 refdcem12 refdcm13]; 
refdcm2-[refdcm21 refdcm22 refdcm23 ]; 
refdcm3-[refdcm31 refdcm32 refdcm33 ]; 
3x3 DCM for reference data 
refdcm-[refdcm] ; refdcm2; refdcm3 |; 


Create reference position and velocity vectors 


refpos=[refx(i); refy(i); refz(i)]; 
refvel-[refxdot(1); refydot(1); refzdot(1)]; 


137 


% Perform transformation from J2000 coordinates to UEN 


refuenp=refdcm*refpos; 
refuenv=refdcm* refvel; 


% Calculate state vector differences 


difp=gpsuenp-refuenp; 
difv-gpsuenv-refuenv; 


vertp(i)=difp(1); 
dtrkp(i)=difp(2); 
xtrkp(i)=difp(3); 


vertv(i)=difv(1 ); 
dtrkv(i)=difv(2); 
xtrkv(i)=difv(3); 


end 

% Calculate average differences 
vertpavg=mean(vertp); 
dtrkpavg=mean(dtrkp); 
xtrkpavg-mean(xtrkp); 
vertvavg=mean(vertv); 
dtrkvavg=mean(dtrkv); 
xtrkvavg=mean(xtrkv); 

% Calculate median differences 
vertpmed=median(vertp); 
dtrkpmed=median(dtrkp); 
xtrkpmed=median(xtrkp); 
vertvmed=median(vertv); 


dtrkvmed=median(dtrkv); 
xtrkvmed=median(xtrkv); 
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% Calculate standard deviations of differences 


vertpstd=std(vertp); 
dtrkpstd=std(dtrkp); 
xtrkpstd=std(xtrkp); 


vertvstd=std(vertv); 
dtrkvstd=std(dtrkv); 
xtrkvstd=std(xtrkv); 


% Plot position component differences 


figure(1) 

plot(gpssec(estr),vertp(estr),'w* ‚gpssec(wstr),vertp(wstr),'wo',... 
gpssec(bdgqc),vertp(bdge),'wx',gpssec(wIrc),vertp(wIrc),'w+',... 
gpssec(kmtc),vertp(kmtc),'w.') 

xlabel(Day 251 Mission Seconds') 

ylabel('Vetical Position Difference (m)') 

grid 

tol=-1; 

legend(ESTR''WSTR''BDQC' 'WLRC''KMTC' tol) 

print 


figure(2) 

plot(gpssec(estr),dtrkp(estr),'w* ‚gpssec(wstr),dtrkp(wstr),'wo',... 
gpssec(bdac),dtrkp(bdac),'wx',gpssec(wirc),dtrkp(wirc),'w+',... 
gpssec(kmtc),dtrkp(kmtc),'w.') 

xlabel('Day 251 Mission Seconds") 

ylabel(‘Downtrack Position Difference (m)') 

grid 

tol=-1; 

legend(ESTR''WSTR''BDQC' 'WLRC''KMTC', tol) 

print 


figure(3) 

plot(gpssec(estr),xtrkp(estr),'w*' gpssec(wstr),xtrkp(wstr),'wo',... 
gpssec(bdqc),xtrkp(bdgc),wx',gpssec(wlrc),xtrkp(wirc),'w+-,... 
gpssec(kmtc),xtrkp(kmtc ),'w.') 

xlabel('Day 251 Mission Seconds") 

ylabel('Crosstrack Position Difference (m)') 

grid 

tol=-1; 

legend('ESTR''WSTR''BDOC''WLRC''KMTC'tol) 

print 
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% Plot velocity component differences 


figure(4) 

plot(gpssec(estr),vertv(estr),'w* ‚gpssec(wstr),vertv(wstr),'wo',... 
gpssec(bdaqc),vertv(bdaqc), wx',gpssec(wlrc),vertv(wlrc),'w-",... 
gpssec(kmtc),vertv(kmtc),'w.') 

xlabel(‘Day 251 Mission Seconds’) 

ylabel('Vetical Velocity Difference (m/s) 

grid 

tol=-1; 

legend(ESTR' 'WSTR''BDQC''WLRC''KMTC' tol) 

print 


figure(5) 

plot(gpssec(estr),dtrkv(estr),'w* ‚gpssec(wstr),dtrkv(wstr),'wo',... 
gpssec(bdqc),dtrkv(bdqc),'wx',gpssec(wlrc),dtrkv(wilrc),'w-".... 
gpssec(kmtc),dtrkv(kmtc),'w.') 

xlabel(‘Day 251 Mission Seconds’) 

ylabel(‘Downtrack Velocity Difference (m/s)') 

grid 

tol=-1; 

legend(ESTR' 'WSTR''BDQC''WLRC''KMTC' tol) 

print 


figure(6) 

plot(gpssec(estr),xtrkv(estr),'w* ‚gpssec(wstr),xtrkv(wstr),'wo',... 
gpssec(bdgc),xtrkv(bdgc),'wx',gpssec(wIrc),xtrkv(wirc),'w+',... 
gpssec(kmtc),xtrkv(kmtc),'w.') 

xlabel('Day 251 Mission Seconds’) 

ylabel('Crosstrack Velocity Difference (m/s)') 

grid 

tol=-1; 

legend(ESTR''WSTR''BDQC''WLRC''KMTC' tol) 

print 
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APPENDIX H. DAY 250 STK PLOTS 
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APPENDIX AA. DAY 259 MATLAB PLOTS 
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